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Abstract 

The Ginibre ensemble of complex random matrices is studied. The com- 
plex valued random variable of second difference of complex energy levels is 
defined. For the N=3 dimensional ensemble are calculated distributions of 
second difference, of real and imaginary parts of second difference, as well as 
of its radius and of its argument (angle) . For the generic N-dimensional Gini- 
bre ensemble an exact analytical formula for second difference's distribution 
is derived. The comparison with real valued random variable of second dif- 
ference of adjacent real valued energy levels for Gaussian orthogonal, unitary, 
and symplectic, ensemble of random matrices as well as for Poisson ensemble 
is provided. 
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I. INTRODUCTION 



Random Matrix theory assumes that the Hamiltonian operator H of a generic quantum 
system is unknown and unknowable [1-3] . The matrix elements of Hamiltonian in given 
basis of Hilbert space are random variables. Their distributions are given by appropriate 
formulae depending on studied Random Matrix ensemble [1-3]. Thy symmetry properties of 
H which is hermitean lead us to Gaussian ensembles of random matrices: orthogonal GOE, 
unitary GUE, symplectic GSE, as well as to circular ensembles: orthogonal COE, unitary 
CUE, and symplectic CSE. The energies E { of quantum systems calculated from diagonal- 
ization of Hamiltonian matrix are random variables with appropriate distributions and 
they exhibit generic classes of level repulsion. It was Wigner who firstly discovered level 
repulsion phenomenon [1-3]. The applications of Random Matrix theory are very broad: 
nuclear physics (slow neutron resonances, highly excited complex nuclei), condensed phase 
physics (fine metallic particles, random Ising model [spin glasses]), quantum chaos (quan- 
tum billiards, quantum dots), disordered mesoscopic systems (transport phenomena). J. 
Ginibre studied very general case of random Hamiltonians. He dropped the assumption 
of hermiticity of Hamiltonians and he considered generic complex valued matrices [1,2,4,5]. 
Thus, i/'s belong to general linear Lie group GL(N, C), where N is dimension and C is 
complex numbers field. Therefore, the energies Zi of quantum system ascribed to Gini- 
bre ensemble are complex valued. This is an extension of Gaussian or circular ensembles. 
J. Ginibre postulated the following joint probability density function of random vector of 
complex eigenvalues Z\, Z N for N x N Hamiltonian matrices [1,2,4,5]: 

N ^ N N 

P{zi, z N ) = H Y[\ Zi - Zj \ 2 -exp(-£ Izj] 2 ). (1) 

j=i 7r 'J- i<j j=l 

We emphasise that Z^s are complex valued random variables, and z^s are complex sample 
points (zi G C). 

One must emphasise here Wigner and Dyson's electrostatic analogy. A Coulomb gas of 
N unit charges moving on complex plane (Gauss' plane) C is considered. The vectors of 
positions of charges are Zj's and potential energy of the system is: 

U{z u ...,z N ) = ~Y J ^\z i -z j \ + ^k 2 |. (2) 

i<j i 

If gas is in thermodynamical equilibrium at temperature T = {{3 = = 2, ks is Boltz- 
mann's constant), then probability density function of vectors of positions is P(zi, zn) 
Eq. (1). Thus, complex energies of quantum system and vectors of positions of charges of 
Coulomb gas are analogous to each other. In view of above analogy one can consider the 
complex spacings A 1 Z i of complex energies of quantum system: 

A 1 Z i — Z i+1 — Z i: i = 1, (AT - 1), (3) 

as vectors of relative positions of electric charges of Coulomb gas. For Ginibre ensemble 
it were calculated the distributions of real valued absolute values of spacings of nearest 
neighbour ordered energies. We complement this by introduction of complex valued second 
differences A 2 Zi of complex energies: 
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A 2 Z t = Z l+2 - 2Z l+1 + Z u i = 1, (N - 2). 



(4) 



The second differences are three energy level magnitudes that enhance our knowledge of 
quantum systems with complex energies. Moreover, A 2 Z i , s can be regarded as vectors 
of relative positions of vectors of relative positions of electric charges. One can observe 
movement of electric charges is Cartesian frame of references or in polar one. Since the two- 
dimensional vectors in Cartesian frame of reference have their projections on co-ordinate 
axes, therefore the real and imaginary parts of A 2 Zj's, namely ReA 2 Zj, ImA 2 Zj, can be 
interpreted as projections of second differences on abscissa and ordinate axes, respectively. 
The radii \A 2 Zi\, and arguments (angles) ArgA 2 Z« of second differences have interpretations 
of polar co-ordinates of A 2 Zi vectors. A 2 Z i , s are analogous to real valued second differences: 

A 2 E, = E i+2 - 2E i+1 + Ei,i = 1, (N — 2), (5) 

of adjacent ordered increasingly real valued energies Ei defined for GOE, GUE, GSE, and 
Poisson ensemble PE (where Poisson ensemble is composed of uncorrelated randomly dis- 
tributed energies) [6-9]. We will calculate the distributions of A 2 Zi, ReA 2 Z;, ImA 2 Zj, 
|A 2 Zj|, ArgA 2 Zj. Finally, we will compare these results with second difference distributions 
for Gaussian ensembles, and Poisson ensemble [1,10-15]. 

II. SECOND DIFFERENCE DISTRIBUTIONS 

We use formula (1) with N = 3 and define the following complex valued random vector 
(Yi,Y 2 ,Y 3 ) and real Aj and imaginary Bj parts: 

Y 1 = A 2 Z l5 Y 2 = Z 2 - Z 3 , Y 3 = Z 3 , Yj = (A„ B 3 ), A 3 = ReY}, B j = ImYjJ = 1, 3. (6) 

The change of the variable formula gives us the result for joint probability density function 
of random vector (Yi, Y 2 , Y 3 ) [16]: 

f(Y 1 ,Y 2 ,Y 3 ){yi,y2,y3) = f(A 1 ,B 1 ,A2,B 2 ,Az,B 3 ){aiMia 2 ,b 2) a 3) b 3 ) = (7) 

! [(ai + a 2 f + (b, + b 2 f] ■ [a 2 + bj] ■ [( ai + 2a 2 ) 2 + (h + 26 2 ) 2 ] • 



12tt 3 

•[exp(-(ai + 2a 2 + a 3 ) 2 - (h + 2b 2 + b 3 f - (a 2 + a 3 ) 2 - (b 2 + b 3 f - a\ - b 2 3 )], 

where yj = (cij,bj) G C are complex random sample points. [In order to obtain Eq. (7) we 
used the complex valued linear map 



(Y 1 ,Y 2 ,Y 3 ) = E(Z 1 ,Z 2 ,Z 3 ),E 



l -2 1 
1 -1 
1 



(8) 



where Jacobian of the inverse map S 1 is equal to unity Jac(S : ) = 1]. We integrate out 
Y 3 : 

/(Yi,y 2 )(yi,y2) = / f(Y U Y 2 ,Y 3 )(yi^y2,y3)dy 3 , (9) 
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and we obtain following marginal probability density function: 

/(Yi,Y- 2 ) (2/1,2/2) = f(A 1 ,B 1 ,A2,B 2 )(a 1 ,b 1 ,a2,b 2 ) = (10) 

[( ai + a 2 ) 2 + (h + b 2 ) 2 ] [a\ + b 2 2 ] {{a, + 2a 2 ) 2 + (61 + 26 2 ) 2 ] • 



36tt 2 



2 2 
• exp[— -a 2 — 2ax(i2 — 2a\ — -fr 2 — 2b\b 2 — 2b 2 ,}. 

Now we calculate the marginal probability density function of second difference: 

f Yl (Vi) = /(A 1)Bl) (ai,6i) = J c f(Y u Y 2 )(yi,y2)dy 2 = ^H( a i + h lf + 24 ] " exp(-^(a 2 + b\)). 

(11) 

Now we derive the distributions of real part Ai and of imaginary part B\ of second 
difference: 

/Ai(ai) = J^/(Ai,Bi)(ai,6i)d6i = 576 ^ («i + 6a 2 + 51) ■ exp(--a 2 ), (12) 



Oijaai = - 

where R is field of real numbers. 

We transform complex valued random variable of second difference Y\ to polar co- 
ordinates' variables i?i,$i: 

^ = 1^1,$! = ^^, (14) 

and we obtain by standard method the following probability density function of random 
vector (R 1: $1): 

Wofa^i) =®(ri)-^rM + 24)- exp(-ir 2 ) (15) 
(the Jacobian of transformation is r±, 9 is Heaviside (step) function [16]). It follows that: 
f Rl (ri) = 0(^)^(^ + 24) -exp(-ir 2 ),/ $1 (00 = e [0,2tt]. (16) 



III. N-DIMENSIONAL GINIBRE ENSEMBLE 

The case of generic N-dimensional Ginibre ensemble is of special physical interest (N > 3). 
We will calculate the distribution of second difference for the ensemble. One commences 
with n-level correlation function: 
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P n (z l , ...,z n ) = P(z!,...,z N )dz n+1 ...dz N = (17) 

jQN-n 

n 

= n- n exp(-J2\zi\ 2 )detD^ n \ (18) 



i=i 

N-l „k 

)('<) • ' > ^ 2 



^ = ejv-i(^),i,j = l,...,n,ejv_i(z) = ( 19 ) 



fc=0 



Ref. [1]. In order to calculate the distribution of complex second difference W\ = A 2 Zi for 
N-dimensional Ginibre ensemble one substitutes n=3 to Eq. (17) and defines random vector 
W = (W U W 2 ,W 3 ) : 



(w 1 ,w 2 ,w 3 ) = n(z 1 ,z 2 ,z 3 ),n 



1 -2 1 
1 
1 



(20) 



The probability density function of random vector W reads: 

P 3 (w 1 ,w 2 ,w 3 ) = P 3 (Q~ 1 (wi, w 2 ,w 3 )) = P 3 (wi + 2w 2 - w 3 ,w 2 ,w 3 ), (21) 
Ref. [16]. Hence, the distribution of second difference is: 

P 3 (u>i) = / P 3 {wi + 2w 2 - w 3 , w 2 , w 3 )dw 2 dw 3 . (22) 
Jc 2 

We combine Eqs (17), (18), (19), (21), (22), and we use Laplace's expansion of determinant 
det £>(">: 

3 3 

P 3 ( Wl ) = or 3 J2(-l) v f exp(- £ l(^~V)i| 2 ) LI e N ^[(Q^w) k ■ {n~ l wy rk ]dw 2 dw 3l (23) 

V ° 2 i=l k=l 

where V is permutation of indices (1,2,3). The only nonzero contribution to Eq. (23) is 
for identity permutation V — id = (1, 2, 3). It results from the fact that other permutations 
produce factors that are periodic functions of arguments Argw 2) Argw 3 of complex num- 
bers w 2 ,w 3 (the integrals over w 2 ,w 3 can be transformed to polar co-ordinates where the 
arguments Argw 2 ,Argw 3 are integrated over [0,2-7r]). Hence, 

3 3 

P 3 ( Wl ) = n- 3 / 2 exp(-^|(fi- 1 H| 2 )ne JV _ 1 [|(fi- 1 «;) fe | 2 ]^2^3, (24) 
Jc ' 2 i=i k=i 

Therefore, considering Eq. (19) one has: 

N-1N-1N-1 i 

P 3 ( Wl ) =n-*J2Y,i: -T-T-l'r .!,,,,(«■>). (25) 
j 1= j 2 =0 j 3 =0 J1-J2-J3- 

3 3 

= L«p(-Ei( n_lw )ii 2 ) n m- i w) k \^dw 2 dw 3 . (26) 
jc i=i fc=i 

One changes variables in Eq. (26) in following way: V 2 = 2W 2 , V 3 = —W 3 , and obtains: 
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= 2' 2 ^ J exp(-|«; 1 + v 2 + t» 3 | 2 - ^\v 2 \ 2 - \v 3 \ 2 )\ Wl + v 2 + v 3 \ 2jl \v 2 \ 2j *\v 3 \ 2j3 dv 2 dv 3 . (27) 

The above double integral can be calculated by extending the exponent by additional terms 
proportional to A« parameters and considering appropriate derivatives: 

Qjl+32+j3 

I nj2J3 ( Wl ) = 2-^ a? . lAi ^ A2ffl3A3 ^i,A 1 ,A 2> A 3 )U <= o, (28) 

F(wi, Ai, A 2 , A 3 ) = / exp[G(tui, u 2 , «3, Ai, A 2 , A 3 )]cfa 2 cfa 3 , (29) 

G(wi,U2, i; 3 , Ai, A 2 , A 3 ) = (Ai - l)|iwi + v 2 + v 3 | 2 + (A 2 - ^)\v 2 \ 2 + (A 3 - l)\v 3 \ 2 . (30) 

Finally, we derive F(wi, Ai, A 2 , A 3 ) by transformation of parametric quadratic form 
G(w±, v 2 , v 3 , Ai, A 2 , A 3 ) to canonical form and integrating over v 2 , v 3 : 

F(w u Ai, A 2 , A 3 ) = A(Ai, A 2 , A 3 ) exp[- J B(A 1 , A 2 , A 3 )|«;i| 2 ], (31) 

where 

" 4(A " A2 ' A3> = (A, + A 2 - |) ■ (a!+ A 3 - |) - (A, - 1)^ (32) 

B(A 1 ,A 2 ,A 3 ) = (A 1 -1).|^^±|. (33) 

Hence, we obtained analytical formula for distribution P 3 (wi) of second difference for N- 
dimensional Ginibre ensemble combining Eqs. (25), (28), (31), (32), (33). What is worth to 
be mentioned is that second difference's distribution is a triple sum of zero-centred Gaussian 
distributions with different widths. The distribution is again function of only modulus |u>i| 
of second difference and it has global maximum at origin. 



IV. COMPARISON 

Finally, in order to compare our results for second difference for Ginibre ensemble with 
previous ones for Gaussian and Poissonian ensembles we must rescale them by division 
by appropriate magnitude. We consider such rescaled dimensionless second differences in 
following way. The mean values of second differences either in real or in complex case are 
zero, hence we cannot divide second differences by the mean values. It follows that we divide 
real valued second differences A 2 ^ for GOE(3) ((3 = 1), for GUE(3) (0 = 2), for GSE(3) 
(f3 — 4), for PE (f3 — 0), by mean spacings < Sp > calculated for those ensembles, and we 
create new dimensionless second differences: 

C,= ^%, (34) 

respectively [6-9]. The probability density functions of (Tg's were calculated for GOE, GUE, 
GSE, and PE [6-9]. Since second difference for Ginibre ensemble is complex valued, then we 
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choose its real part A\ for comparison with Cps. One divides A\ by analogue of < Sp >'s, 
which is mean value < R ± > of radius i?i Eq. (16) of A 2 Z X : 



<i? 1 >=y o nfn^dn = — Vto. (35) 



Hence: 



• Y ' - T^P < 36 » 

is rescaled dimensionless Ai. The probability distributions of Cp's and of Xi depend on the 
same real variable x which is equal to , < ° 1 i> , respectively (ei is value of energy i^). 

The second differences' distributions for Gaussian, Poisson, and Ginibre ensembles as- 
sume global maxima at origin and that they are unimodular. Firstly, it extends the theorem 
of level homogenisation to Ginibre ensemble [6-9]. We can formulate the following law: En- 
ergy levels for Gaussian ensembles, for Poisson ensemble, and for Ginibre ensemble tend to 
be homogeneously distributed. The second differences' distributions assume global maxima 
at origin no matter whether second differences are real or complex. From Coulomb gas' 
point of view it is easier to be understood. The unit charges behave in such a way that the 
vectors of relative positions of vectors of relative positions of charges statistically tend to be 
zero. It could be called stabilisation of structure of system of electric charges. The above 
results can be extended to study of higher differences' distributions for Ginibre ensemble. 
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